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Abstract 

The self-consistent theory of the fermion condensation, a specific phase transition 
which results in a rearrangement of the single particle degrees of freedom in strongly 
correlated Fermi systems is developed. Beyond the phase transition point, the single 
particle spectra are shown to be flat. The interplay between the flattening and the 
damping of the single particle spectra at T ^ is investigated. The width 7(e) of 
the single particle states is found to grow up linearly with e over a wide range of 
energy as in a marginal Fermi liquid. Our results gain insight into the success of the 
phenomenological theory of the normal states of high-temperature superconductors by 
Varma et al. 

1. Introduction 

In the article |]1| within the Landau-Migdal (LM) quasiparticle pattern, a new phase 
transition in strongly correlated Fermi systems, the so called fermion condensation, has 
been uncovered. Its salient feature is in the appearance of the fermion condensate (FC), a 
group of degenerate states whose energy at T = coincides with the chemical potential /i. 
Owing to the degeneracy, the quasiparticle occupation numbers n(p,T = 0) are no longer 1 
or 0. They are determined by a variational condition [|l|, 0, |^ 

(5Eo(n)/5n(p) =/i, pen, (1) 

where EQ{n{p)) stands for the ground state energy, while fl denotes the FC region whose 
boundaries are determined by eq. (P itself. According to Landau, the l.h.s. of (0) is nothing 
but the quasiparticle energy e{p) and therefore this equation implies smearing the Fermi 
surface — its metamorphosis it into a volume in the three-dimensional systems or into a 
surface in the two-dimensional ones. In anisotropic systems, say, in crystals, the FC may 
occupy only part of the old Fermi surface, e.g. a region, adjacent to the van Hove singularities 

Bi- 

At T > the quasiparticle distribution n(p, T) is found with the help of the usual Fermi- 
liquid-theory variational equation 6F/6n{p) = |^ where F = Eq — fiN — TS is the free 
energy. The solution is written in the standard Fermi-Dirac form 



n(p,T) 



1 + exp 



T 



(2) 



but the energy ^(p) = ^(p) — fi depends on the distribution n itself. For both equations 
(E)) (i) to be consistent, the spectrum ^(p) should exhibit a strong T-dependence in the FC 
region 

eMp,T)^Tlnl^^, pen, (3) 
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where no(p) is the solution of eq. (|I]). We see that the plateau in the spectrum ^(p) existing 
at T = is converted into an inclined plane with the slope proportional to T. Such a 
significant liability of single particle characteristics with respect to T is a calling card of 
fermion condensation. 

The fiattening of the single particle spectra ^ (p) has been experimentally observed in the 
ARPES study of electron systems of some high-Tc superconductors p, 0, H but observable 
single particle peaks have much larger widths 7 than those calculated employing Fermi liquid 
theory. As a result, the applicability of the quasiparticle formalism to such systems has been 
put in question and alternative explanations have been proposed |10 |. 



As known in Landau theory, the single particle spectrum ^{p) is described by the only 
parameter — the effective mass M*, characterizing the slope {d^{p)/dp)p ~ (M*)~^ of the 
spectrum at the Fermi surface. Usually the effective mass is treated as a T-independent 
quantity. But, as seen from @, beyond the FC phase transition point, the slope d^pc/dp 
is just proportional to T, i.e. for the FC states M* ~ T~^. This peculiarity of the fermion 
condensation alters the Landau theory predictions even if the quasiparticle picture survives 
in the systems with the FC. However, the ratio •ypciT) /C,pc{p,T) calculated within the 



Landau approach ||TT| |T2|, |T3[ does not vanish at T ^ in contrast to fundamentals of the 
Fermi liquid theory. In other words, at finite T, the theory of the fermion condensation 
based on the LM quasiparticle concept suffers from inconsistency. 

On the other hand, at T = 0, the quasiparticle picture is known to hold as long as 
the critical index u characterizing behavior of the imaginary part of the mass operator 
Im E(p, e) ~ e'^, exceeds 1 (here e is measured from the chemical potential fi). This conclusion 
stems from the formula for the renormalization constant z given by 



1 - Re 



fdJ:{pp,e] 



(4) 



V de 

Straightforward calculations with taking into account the well known dispersion relation 



14, 15 



demonstrate that the factor z vanishes provided u < 1. Thus if z/ < 1 the conventional 
quasiparticle picture is destroyed. We shall see that this just occurs beyond the point of the 
fermion condensation. 

The interplay between damping and flattening beyond this point is quite subtle. Indeed, 
let the solution with the FC, i.e. with the plateau in ^(p) at T = 0, exist neglecting damping 
effects. When damping effects get involved the plateau in ^ (p) should somehow change. If the 
flat portion of ^ (p) were completely destroyed then the spectrum ^ (p) could be characterized 
by the single constant M* that results in reviving Fermi liquid theory and the FC solution, 
as well. The way out is that the flat portion in the single particle spectrum survives but 
flattening and damping should somehow counterbalance each other. 

In this article, we develop the self-consistent theory of the fermion condensation proceed- 
ing from the standard many-body-theory approach and bypassing the variational condition 
(P which is not adapted to take damping effects into account. The flrst step is to derive 
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a general relation between the spectrum ^(p) and the distribution n(p) assuming the LM 
quasiparticle pattern to be valid. The second step is to verify that this relation considered 
as an equation has a class of peculiar solutions for which ^{n{p)) = in a whole region 
fl of the momentum space rather than at several isolated points. These are just the FC 
solutions. However, as we shall see, the LM quasiparticle picture is nothing more than a 
reasonable approximation for the treatment of the systems with the FC. These systems do 
possess features inherent in marginal Fermi liquids for which the quasiparticle weight in 
the single particle state adjacent to the Fermi surface approaches zero. The next step is 
to extend the Landau formalism to marginal Fermi liquids including the derivation of the 
relation between the single particle spectrum and the distribution function. Then one needs 
to calculate irregular contributions to the single particle mass operator beyond the FC phase 
transition point and evaluate how strong the plateau in the spectrum of the single particle 
excitations is inclined due to the damping. 

2. The fermion condensation in the quasiparticle approximation 

We already mentioned above that the systems with the FC possess the properties inherent 
in marginal Fermi liquids. In particular, the single particle mass operator E(p,£:) contains 
irregular terms precluding expansion in a Taylor series in the energy e. However, these terms 
rapidly die out if e exceeds a characteristic FC energy epc- Therefore the quasiparticle 
picture seems to be a reasonable first approximation for studying the energy dependence of 
quantities of our interest. 

In the system with the FC, the quasiparticle Green function is found by retaining 
leading terms in the expansion of the inverse Green function in the limit £ — > 0. This 
expansion reads as 

G-\p,e^ 0)^e-e,- S(P,£=0) - e (^%^)^, (6) 

where = /2M — fi. Here we neglected the variation of the derivative (9E(pc, £)/de)Q in 
the FC region and calculated this derivative at a fixed point Pc G ^• 
The pole term has the accustomed form 

However, the quasiparticle spectrum ^(p) near the Fermi surface is given by a different 
formula 

e(p) = ^(e; + s(p,£=o)), (7) 

than in the Landau theory of normal Fermi liquid since the momentum expansion with 
retaining terms of the type {p'^—pp){dIl/dp'^)F makes no sense in the system with the FC. 
The quasiparticle momentum distribution is evaluated as usual 

n(e(p)) = Im/G^(p,£)^, (8) 
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yielding the step function 6'(^(p)). 

Another relation between ^ and n is found proceeding from the Landau-Pitaevskii identity 



Em 

Here ?7 stands for the block of diagrams irreducible in the p-h channel. 

In doing the renormalization of eq. (|^) a little modification of the standard decomposition 
l^, IT5[ of the product GG into a regular part B and a singular one A is necessary. It refers 



to the singular term A. Now it is defined as follows 

G^(p,£)G'?(p+q,£+o;)— . (10) 

Then 

f dG-\p,6) de _ 2 r gG-i(p,e=0) & _ an(p) 

i ^^^P' 2^ - ^ i ^ ^P' 9^:^ ^ ^p' "^2^ - 

In obtaining this result, we have used relation 

dpa dpa 



As a result, after conventional transformations, we arrive at the well known relation |TJ, p!5 

^^(P.^) Pa . [ an(pi) , i^^. 

The phenomenological quantity JF ~ (1 — BU)^^U, being a;- independent, plays the role 
of an effective interaction potential between particles in medium. Assuming this function to 
be known including possible dependence of on the distribution n{^) itself, eq. (^TJ) can be 
integrated and conveniently rewritten as 

2 

^^P^ ^2M^I ^^P' P'' ""^^^P^^^ ^^^^ 

the effective Hartree-like term H{p, pi) being related to the effective interaction JF by the 
formula 

on 

As mentioned above, the momentum distribution function n{^) is simply 6{$,) so that eq. (0) 
can be treated as an equation for finding the spectrum ^(p). Usually and, hence, H are 
nearly momentum independent functions. Then eq. ([12D has the single solution ^ (p) for any 
p, with d^/dp > 0, i.e. there is one-to-one correspondence between ^ and p. 

The situation drastically changes in strongly correlated systems where momentum de- 
pendent components of the effective interaction can be strong. The case of homogeneous 
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matter where particle interact by means of limited long range effective forces, i.e. H{p, pi) = 
VS{p — Pi), analyzed in 0] is especially impressive. In this case, relation between ^ and the 
distribution n becomes algebraic: 

^ = C + V^(0- (13) 
This formula considered as an equation for finding C, {p) is rewritten in the form 

L{^,p)^^-^=n{0, (14) 

where n{^) = 6{^). This equation has a simple graphical solution (see Fig. 1) that furnishes 
a transparent explanation of how the FC solution of eq. (^2|) arises in more complicated 
cases. Indeed, the l.h.s. L{^,p) of eq. (|14D depending explicitly on the momentum p as a 
parameter is depicted by a set of straight lines with the slope ~ 1/V. On the other hand, 
the r.h.s. n{^), made up of two horizontal lines with the kink at ,^ = 0, is independent of p. 
Crossing points of the set L{C,,p) with the horizontal pieces of n{^) are of no interest while 
in the entire interval Pi < p < pf satisfying the inequality 

- r < < 0, (15) 

the straight lines cross the vertical line ^ = 0. Thus all particles with the momenta p varying 
in this interval turn out to have the same energy ^(p) = 0, i.e. they belong to the FC region. 
It is worth noting that the FC solution is stable provided V > 0. 

The graphical procedure is easily generalized to the case of nonzero temperatures where 
the distribution function n(p) is given by eq. (2). This model is also applied to anisotropic 
systems 0] where the energy becomes angle-dependent that results in the anisotropy of 
the FC occupation. So, in a two-dimensional tight-binding model with the single particle 
spectrum given by 

= —eo{cospx + cos Py + t cospx cos Py) — /i, (16) 

and t > 0, the condition (plSf ) of the model for particles becoming trapped into the FC reads 
as 

u V 

< cos -|- COSPy + t COS COSPy H < — . (17) 

It is worth noting that often components of the effective interaction between electrons in 
solids have quite narrow peaks. So, the spin-spin component considered in the first quadrant 
of the quadratic two-dimensional Brillouin zone has the peak close to the point Q = (tt, tt) 
associated with antiferromagnetism WTj. The resonance part of this interaction can be 



crudely written in the form 

5(pi, P2, q) = -A aia2 0(pi) 0(p2) S{q-Q). (18) 

This is an exchange spin fluctuation term. Inserting it into the RPA formula for the mass 
operator E an essentially momentum dependent contribution 

S(p) = ^</.(p)</)(p-Q)n(p-Q) (19) 
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arises (see Fig. 2). Assuming other contributions to be reduced to a renormalization of the 
chemical potential a system of two equations 

ei = e°(p) + ^0(p)0(p-Q)M6), 

6 = e°(p-Q) + ^0(p)0(p-Q)n(ei) (20) 

is obtained. We notice that the graphical procedure of solution requires operating in the 
three-dimensional space. Analyzing the landscape ^(p) we see that the interaction term 
does not affect the region close to the center (0,0). At the same time, its involving gives 
rise to an elevation of the area adjacent to the van Hove points that eventually results in 
the appearance of the FC in this region since otherwise energy of the occupied states would 
exceed that of unoccupied ones. Results of numerical calculations are given in Fig.3, where 
the shadded areas correspond to the regions of the FC. 

The models with limited long range forces are not exceptional. Similar results are ob- 
tained, e.g., in a model with finite range forces given by the Yukawa formula in the momentum 
space: H{pi — P2) = g exp{—P\pi — P2I) p|, H. Now the l.h.s. of eq. (|14D is recast as 

L(e,p) = (-A + /?2)^ = n(0. (21) 

In this case, the graphical representation of solution is complicated by the presence of the 
term A^. Fortunately, in the FC region, this term vanishes and we come back to eq. ([T^). 
Some complifications emerge because of the existence of discontinuities in the term at the 
boundaries of the FC region giving rise to jumps in the distribution n{p) at these points. 
The range of the jumps, which drops with increasing the coupling constant g, at relatively 
small g may attain 1 and then the FC solution ceases to exist. As a result, requirements 
for the appearance of the FC solution become more severe: the constant g should exceed a 
nonzero positive critical value gc- 



To make the next step, let us recast eq. (|12D as follows 

J H-\p, pi) e(pi) dn - Xo(p) = n(0, (22) 

where Xq = (if~^^p). The model ( pT]) is recovered provided H^^{p,pi) = g~^{[3'^ — A) 

5(p — pi). For a certain class of the momentum dependent operators this expansion 
can be extended to 

^"'(P, Pi) ^(Pi) = (^0 - ^2A + A4A2 . . .) 5{p - pi). (23) 

Upon substituting this equation into eq. (0) one finds 

(Ao - A2A + A4A2 . . .)(e - H) = niO- (24) 

Again all the terms containing the Laplace operators applied to ^ vanish in the FC region 
where ^ = 0, and the FC solutions ceases to exist provided the jumps of n at the boundaries 
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of the FC region attain the unity. Thus we infer that for a quite wide class of the effective 
single particle Hamiltonian if, the FC solutions ^ = exists. The extension of this result to 
the general form of the operator i7(p,pi) is straightforward. We seek the FC solutions in 
an equation 

= + 1 H{p, pi; n) n(e(pi)) dn, pen, (25) 

obtained from eq. (|T2p setting there ^ = 0. 

It is the usual integral Fredholm equation of the first kind. As known, numerical solving 
such equations encounters difficulties. Specialized methods have been developed to overcome 
these difficulties which are often called inverse problems (see e.g. ||20[). We are going to 
analyze them in a separate paper. 

Thus we see that the FC solutions can be obtained proceeding from the Landau-Pitaevskii 
relation between the single particle spectrum ^(p) and the quasiparticle momentum distri- 
bution n(p), without addressing to the variational condition (|^). What is remarkable is that 
the Migdal jump in the distribution function n at the point ^ = 0, being a building block of 
Fermi liquid theory in systems with velocity independent forces, becomes a principal factor 
for the appearance of the FC solutions in systems where such a dependence turns out to be 
sufficiently strong. 

To conclude, having at hand the FC solution ^(p G fi) = 0, one can recover the variational 
condition (|l]) suggesting e{p;n) = 6Eo{n)/6n{p) and remembering that ^(p) = £:(p) — /i. 
In principle, the functional Eo{n) can be found from this variational equation provided the 
dependence of e(p, n) on the distribution n is supposed to be known. 

3. The extension of the Landau formalism to marginal Fermi liquid 

If the systems with the FC were subject to the Landau theory then the development of 
the self-consistent theory of the fermion condensation would be completed. However, this 
is not the case. As we shall see, the mass operator S(p, e) in these systems contains the 
irregular logarithmic term ~ eln|£| inherent in marginal Fermi liquid. This term makes it 
impossible expand ^{e) in a Taylor series in energy. Although the pole in the single particle 
Green function G survives but the residue turns out to be logarithmically small vanishing 
at the Fermi surface (see Fig. 4). 

In spite of this fact, there is good reason to believe that underlying ideas of the Landau 
theory remain robust including that of the decomposition of the single particle Green function 
G into a sum G = G"* + ZG^ with a singular part G^ (we shall call it the pseudoparticle 
propagator) containing the logarithmic components of the mass operator. In homogeneous 
matter, we write in a simple form 

G' {P,e) = , , ^ ' . , , . (26) 

e — a{e) — s{p) — tK{€) 

The connection between the quantities cr and n stems from the dispersion relation (13) which 



implies ^1 



nie 0) = —KS, ijie — > 0) = —ne In — , (27) 

TT El 
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where k and El ^ B^p = p'jp/2M are numerical constants. 

In the Landau theory, the quasiparticle distribution n(^) is the universal function with 
the kink at ^ = 0. In marginal liquids, the pseudoparticle distribution function n{^) ceases 
to be universal. Analogously to eq. (8), it is evaluated from relation 



1 °° 

n(0 = - / lmG'ip,e)de 

IT J 



1 1 



'y{e)de 



'y{e)de 



2 



:28) 



where we introduced notations ^ = Zs,'j = Zk and ctl^s) = ;^7£:ln^. The constant Sl 
in (^) can be chosen in such a way as to ensure the equality between the particle and the 
"pseudoparticle" numbers, i.e. 



Tr J G'{p,e) dr 



de 
2Td' 



(29) 



We notice that if the ratio 7(^)/^ is relatively small then the integral (^) is close to 1, 
provided ^ < 0, or to in case > 0. On the other hand, if this ratio is of order of unity 
(marginal regime) then there n(^) essentially deviates from ^(^). As a result, the strength of 
the FC states turns out to be distributed over all this interval and it is its upper boundary 
Epc that should be compared to the characteristic FC energy. 

The renormalization of eq. (P) in the marginal region is implemented along the same 
lines as before with the natural redefinition of the singular part A which now contains the 
product C^C^. Furthermore in contrast to the quasiparticle approximation, ^(p) does not 
coincide with the spectrum of the single particle excitations determined by the location of 
the pole of G'*(p, e) although both these quantities are related to each other unambiguously. 
After some algebra one finds 



dpa 



Pa 

M 



+ 



J -^(p,Pi 



^ dn{pi 
dpia 



-dri. 



(30) 



In obtaining this result, we ignored the possible alteration of the regular part B due to 
corrections related to singular terms in the mass operator. 
Again the integration of this equation yields 



2 

C(P) = ^ + / ^(P' Pi) ^(Pl) 



(31) 



However, in marginal Fermi liquid, the dependence of the distribution n(^) on ^ is intricate 
and for finding the flat portion of the spectrum, one needs to solve eqs. ( pTj ) and ( pH]) self- 
consistently. The result crucially depends on the magnitude of the marginal term in the 
mass operator E. 
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4. The evaluation of the irregular part of E beyond the phase transition point 



The purpose of this chapter is to calculate the imaginary part of the mass operator S(p, e) 
in the presence of a flat portion in the single particle spectrum. In doing so, the temperature 
T is supposed to exceed a temperature Tc of a phase transition, say. Cooper pairing or crys- 
tallization, that results in lifting the degeneracy of ^ (p) implied by the variational condition 
(|1]). In this article, we investigate behavior of the mass operator S(p,£:) in the limit T -C e 
and therefore set T = 0. Without loss of generality, we restict ourselves to the decay of a 
particle implying e > in all the formulas. Different contributions to lmS(p,£: > 0) fall 
into two categories: i) the particle decay into a particle and a particle-hole pair and ii) the 
decay into a particle and a collective state. The evaluation of the last contribution requires 
special calculations of the collective spectra in the systems with the FC. This problem is 
beyond the scope of this article. It will be studied later. With these restrictions, ImS(p,e) 
is given by []T1|, |TB| 



Im E(p, e>0)= JJ drdr-s j j ^(P' Pi' ^i' Pa, ^2, Ps, ^s) 



X ImG'(pi,£i)ImG(p2,e2)ImG(p3,£3)- (32) 

The 4-momenta p,p2 refer to incoming particles while pi,P3, outgoing ones: pa = P+P2— Pi; 

= e + 82 — El = 62 + uj. The limits of integration over £3 are imposed by the restrictions: 
i) £3 > since the 4-vector ^3 corresponds to the particle and ii) E3 < u since E2 < because 
the 4-vector p2 corresponds to the hole. In the three-dimensional system dr = d^q/ (27r)^. But 
most of experimental data on properties of strongly correlated Fermi systems are obtained in 
studies of anisotropic quasi-two-dimensional crystals where dr = (i^g/(27r)^. The transition 
probability W is related to the square of the effective interaction amplitude between particles 
averaged over spin variables |T8[ 

Ty = l|rTtP + ^|rnp. (33) 

Our prime goal is to establish the energy dependence of the mass operator S(p, e) beyond the 
FC phase transition point and further numerical factors will be omitted. In this article, we 
also neglect terms containing C since their contributions are reduced to a renormalization 
of input parameters. 

Usually, the interaction F varies slowly versus its arguments. As a result, the integration 



is essentially facilitated [|T^ . But it is not always true. We shall see that in the systems with 
the FC, behavior of F as a function of u drastically changes that has a dramatic impact on 
the energy dependence of the mass operator S. At the same time, dealing with the damping 
effects the momentum dependence of F is of no crucial importance. Therefore to avoid 
unjustifled compliflcations we restrict ourselves to the momentum-independent zero Landau 
harmonic of the interaction amplitude. Then F^^ = while T^i is expressed in terms of the 
scalar component F of the interaction amplitude F|| = 2F. This suggestion being consistent 
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with our choice ( p6[ ) of the form of furnishes the opportunity to integrate over ps and 
rewrite (13^) in the form 



e 



^ ' J I 27r[£-a;-a,.(p-q,£-^)-e(p-q)]2 + 72(p-q,^-^) ^ ^ 

where q = p — Pi- 

The interaction F is related to the effective interaction potential JF by the well known 
Landau equation 

r(q, uj) = ^(q) + ^(q) A(q, uj) r(q, u;) = f^^j -, (35) 

1 -^(q)^(q,^) 

developed in the particle- hole (p-h) channel. As usual, the factor is already included into 
the definition of the effective interaction T. Then the propagator A{c[,uj) given by 

de 



A(q, uj) = 2[ G^(p, e) G^(p-q, e-uj) dr (36) 
J 2m 



is conveniently rewritten |14 



as 



A{q, Lu) = NniO) (a(q, tu) + ib{q, 



ijj 



with A^„(0) = "ppM* j'n'^ being the density of states of the normal Fermi liquid and M*, the 
effective mass of the normal quasiparticle, and 

1 /" /" de 
Kq,^^) = ^^p^yy ImG"(p,e)ImG'"(p-q,e-^)rfr— , 

a(q,c^) = ^^//lmG^(p,e)Re[G^(p-q,£-u;) + G^(p-q,e+u;)]rfr^. (37) 

— oo 

In ordinary Fermi systems, the magnitude of |yl|/A^„(0) is of order of 1. 

In the strongly correlated Fermi systems, such as liquid He-3, dense neutron matter or 
electrons in high-temperature superconductors the amplitude T is repulsive [T > 0) and 
strong so that, as a rule, the term \J-'A\ in the denominator of (pSf ) significantly exceeds the 



unity. Omitting 1 in eq. (B^), we are left with [IT9 



In this limit, the interaction amplitude F is completely independent of characteristics of the 
effective interaction potential JF. It holds over a wide range of frequencies being violated 
only close to lines ujs{q) corresponding to the collective excitations of the system. 

Irregular contributions to S may also emerge in the particle-particle (p-p) channel, that 
is important, e.g. in superfluid systems. Here we study the systems for which the critical 



10 



temperature Tc of the superfluid phase transition is very low (e.g. Sr2Ru04 with Tc ^ 1 K) 
and consider the case T > T^, at which the p-p correlations make no difference and therefore 
below we omit these contributions. 

In what follows we consider a crystal with a cubic or square lattice assuming the FC to 
be positioned in the vicinity of the van Hove points 0, 0. We also set that the FC density 
Pc = VP is rather small: rj < 1. The total damping 7 from (0) is decomposed into a sum 

7(P, = 7o(£) + 7i(P, + 72(P, £)■ (39) 

The term 70 does not contain the final FC Green functions. Contributions containing the 
single final FC Green function are included into the term 71, those with two FC ones, into 
the last term 72. It should be indicated that for the square or cubic lattice the decay of the 
FC particle into the final state with no less than two FC states unambiguously results in the 
appearance all three FC final states. 

Dealing with the momentum dependence of 7(p), two regions can be distinguished: i) the 
FC region Q where the single particle spectrum e(p) is flat and ii) the ordinary region where 
the gradient d^/dp^ has the usual order of value ~ pp/Mq where Mq being independent 
of T is the effective mass of the ordinary quasiparticle. Further we shall neglect a possible 
variations of the relevant quantities inside these regions and use the subscript c for the FC 
region and the subscript n, for the "normal" one, e.g. 7(p G i7) = 7c or 7(p ^ Q) = 7„. We 
also use the notation Foo for the set of diagrams of the interaction amplitude without the 
FC states and the notation An{ci,uj), for the particle-hole propagator of the normal Fermi 
liquid. 

In many respects, the evaluation of the quantities To(){q,u) and An{q,uj) at g ~ pi? 
is equivalent to that in normal homogeneous Fermi liquid. The momentum integration in 
Feynman diagrams is carried out with the help of the standard formula 

dp^dpydp^ = dSdpy = .j^f'} | , (40) 

\di/dpy\ 

where v indicates the normal to the Fermi surface direction. 

Overwhelming contributions to 70 come from the momentum region quite far from the van 
Hove points. Here the components 60 and of the propagator An{q,uj) and the function 
|roo(5', t^)P are known to vary slowly enough that allows us to replace them by averaged 
values: 



Upon substituting (^) into eq. (0) the integrals cease to depend on p at all, and we arrive 
at the common result 

7o(^) = (42) 
Pf 

where 70 is a dimensionless positive constant. 

When evaluating 7i(p,£), two ordinary Green functions entering eq. (|32D are combined 
into the particle-hole propagator if the FC particle is generated, or into particle-particle 



11 



one when the FC particle is annihilated. As before, we concentrate on the analysis of the 
contribution 7^^. Then the contribution 7i„ is given by (see Fig. 5) 

e 

7i„(e) ~ A^„(0) / / dTid£ilmGc{vi,^i) \Tm{Tp,V>i,^-^i)? bn{e-ei), (43) 



n 



where 



ImG'c(pi,ei) = ^ ^ -. (44) 

[si - ac(pi, £1) - ^c(pi)J^ + 7c (Pi, ^1) 

The region f2 of the momentum integration in this formula consists of "patches" of the FC 
adjacent to the van Hove points. 

It can be verified that the interaction amplitude Fqi with the single (final) FC state has 
the same order of value as Foo ~ 1/A^n(0). Then eq. (|^) becomes 

M* f f 

In obtaining this result the relation ( ^T]) for the averaged value 60 was used. 

The integral ( ^Sj ) considered versus energy e, the upper hmit of integration, is, apparently, 
a monotonic function. If the upper limit were extended to infinity, the integral value would 
be equal to (1 — Tio(p)) provided ^(p) > 0. In the opposite case ^(p) < 0, it is close to zero 
(recall, we consider the particle decay). Integrating over momenta, one finds 

M* r 

The evaluation of this integral is facilitated if the energy e exceeds the characteristic FC 
energy epc- In this case, in contrast to the density of the states No{e), which varies slowly 
with e, the density of the FC states Nc{e) does extinct at e > epc- Therefore it can be 
approximated by the 5-function: Nc{e) = PcS{e). As a result, at e > epc the FC contribution 
to 7„ turns out to be linear function of energy rather than the quadratic one as the ordinary 
term 79: 

7i„(£ > epc) - -lin V crin(e > Bfc) - £ In — , (47) 

71 El 

where 7i„ is a constant of order unity. 

As for the corresponding contribution 7ic(£^) to the FC damping (see Fig. 6a), the equation 
for its evaluation has the form 

e 

7ic(£) ~ A/'„(0) J I rfrirf£iImGc(pi,ei) |r2o(p,Pi,£-ei)P&n(e-ei), (48) 
n 

where stands for the amplitude of the transition of the normal p-h pair into the FC pair. 

The energy dependence of the quantities in this integral is more important than the 
momentum one and therefore, in the first approximation, the amplitude r2o(p, Pi, £—Ei) can 
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be replaced by the averaged over the FC region value r2o(£^— Then after implementation 
of the momentum integration one finds 

£ 

7i,(£) ~ iV„(0) j N,{e{)\T2^{e-e^)\Hn{e-e^)de^. (49) 



Two final FC states also come into play when the contribution 72n(^) to the damping 
7(e) is considered. Again we concentrate on the p-h contribution (see Fig. 6b) given by 

e 

l2n{£) ^ Nn{Q) j j rfrirf£:iImG„(pi,ei) |r2o(p,Pi,e-£:i)|^&c(p-Pi,e-£i)- (50) 
c 



We introduced a special notation Vt' for the integration region in (|50|) , determined by the 
requirement for the difference p — pi to have the value appropriate for the pair of the FC 
states. Neglecting the momentum dependence of the amplitude and the propagator 62 
this equation is recast to 

72„(£) ~ iV„'(0) j \V2^{e-e{)\H,{e-e{)de^. (51) 


At last, deahng with the contribution 72c(£), all four states belong to the FC (see Fig. 6c), 
so that 



'l2c{e) ~iVn(0) j j c?rid£:iImGc(pi,ei) |r22(p,Pi,e-£:i)P&c(p-Pi,£-£i)- (52) 
n 

This integral is reduced to 

e 

l2c{e) ~ Ar„(0) / N,{e{) \V22{e-e{)\H,{e-e{)dei. (53) 



The averaged over the FC region propagator hc{oj) entering this formula is defined in 
terms of 



1 

lc{e) 7c(e-^^) de 



n 



dr 



(54) 

Primarily, we study the case e > Spc- In this situation, the integrand in eq. ( |5^ ) is made 
up of two separate peaks. The first is positioned at e ~ epc and contributes only if .Cc(pi) > 
(i.e. if Tio(pi) < 1/2) independently of the sign of the other term ic{vi)- The second is located 
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at {u — e) ~ epc and contributes provided ^c(p— q) < (i.e. no(p— q) > 1/2) independently 
of ^c(Pi)- Both these terms supplement each other and we arrive at 



bc{uj > Efc) 



VP 



7c (t^) 



7c (e) de 



7c (5 -t^) de 



[{e-uj- ac{e-uj)f + -il{e-uj) 



+ 



(55) 



Again we can use the fact that the integrands on the r.h.s. of this equation are 5-like functions. 
Then after simple integration we are left with 



where 



(56) 



(57) 



To correctly estimate the damping one also needs to know the behavior of the real part Qc 
of the propagator averaged over the FC region. It can be traced from the formula 



iVn(O) 
1 



fi —00 



ImG'*(pi,e) Re G"'(pi-q, e-cu) + G"'(pi+q, e+cu) 
7c (e) de 



dr de 



dr 



X 



Q —00 

e—uj—o-r(e—u!) 



(Trie 



+ 



e+LU—aAe+LU] 



\e-uj-acie-uj)-^c{Pi-(lW+'yci^-^) [e+uj-(Tc{e+cv)-^c{Pi+(lW+'yci^+^ 
After simple algebra, the following estimate for the averaged value ac{uj) is obtained 



\ac{uj)\ 



VPf ^fc 



M* [cu - (Tc(^)]2 + 7c^(^) 



• (58) 



(59) 



Comparing jad to \bc\ and bearing in mind that the value of epc depends by itself on 7c we 
infer that, in the first approximation, the Cc contribution can be neglected. 

It is worth noting that the integrands of ( ^9]) and (|^) are peaked at ei ^ epc due to 
the (5-like behavior of ImG'c(pi,ei) and rather smooth behavior of the products 6n|r2oP and 
^c|r22| • To verify this assertion consider the set of equations for the amplitudes T22, ^02- In 
the symbolic form, it reads as 



ro2 



-^22 + ^22^cr22 + ^20^11^02 
^02 + ^ 32^cr22 + ^ 30^nr02- 



(60) 
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This system is easily solved 



-^20 

:i-J-00^n)(l-J^22Ac) - J^loA^. 



Y _ -^22(1— -^OO^n) + ^20^n ^^^-^^ 

From the last formula we can infer that in the strong coupling limit the value |Acr22| is 
of order of 1. In what follows we employ an interpolation formula 

iVn(0)|r22(^)r~^^^. (62) 

which is correct for a moderate interaction and in the strong coupling limit, as well. As 
for the ratio |r2o/r22| ~ ^2o/-^22, its value is, probably, considerably less than 1. Upon 
substituting these results into eqs. (^Qf ) and (|53|) we are led to the following equations 



As a result, the equation for finding 7c(£) takes the form 

7o(.) = -0.^ - C , + A „p% J^^^^y (64) 



Upon substituting here the formula (|56| ) for the propagator be and simple transformations 
one obtains 

where 5p = rjpp/M* while the constants A, B, C are of order of 1. At £ > dp, the magnitude 
of the propagator be is less than 1 and therefore the second term in the brackets of eq. (|65|) 
is relatively small. This implies 

7c(£ > 5f) ~ -C7]e, (66) 

where the constant C is of order of 1. We infer that at these energies the magnitude of the 
marginal terms in the mass operator is quite small. It should be indicated that, according 
to (pSD, the presence of the marginal term in S(£ — 0), no matter how small it is, alters the 
Landau distribution n(^) giving rise to the an inclination of the plateau in the single particle 
spectrum proportional to 7] if this behavior persists in the limit e 

However, this smallness cannot always hold as e together with f{e) go down to zero, 
otherwise the sum of the terms in the brackets on the l.h.s. of eq. ( |65D would change its sign. 
This circumstance makes the function 7(5) rapidly enhance so as to both the terms in the 
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brackets almost cancel each other approaching e ~ Sp- Indeed, let us seek the solution of a 
simplified equation 

1 

It can be easily verified that no solution of this equation exists at e exceeding a critical value 
^cr = Spy/A and 7c(£cr) =0. As e goes down, 7c(e) grows up so that the ratio = 7c(£)/£ 
attains in the region e ^ 6f values of order of the unity in contrast to perturbation theory 
estimate — t] <^ 1. The alteration of behavior of the damping from the perturbative 
regime to a "strong coupling limit" 

-fc{e) ^e, e < Sp, (67) 

occurs in a quite narrow region of energy that is confirmed in numerical calculations. 

The result VciSp) ~ 1 is crucial for finding out behavior of 7c(£: 0) where a scaling 
approach can be employed. We assume that 7c (e ^ 0) ~ e:*^. For the evaluation of the 
critical index u, we calculate the propagator be proceeding from an assumption that the 
damping term exceeds all the rest ones in lmGc{p,e). Then 

bc{uj) r-^rj- ^ 



This function diverges at uj —>■ provided u > 1/2 while the product 



\T22{e)\'b,{e) 



r]uj 



Then the folding with lmGc{p,e) ~ l/7c leads to the strong-coupling-limit equation for 7c 
which reads as 

7c(£ ^ 0) ^ / 7c(^) — • (68) 



Significantly that rj does not enter this equation. This allows us to find the correct solution 
7(e) ~ £ which simultaneously meets eq. (|68D and the matching condition 

rddp) ^ 1. 

Thus the marginal behavior of the mass operator S(£) in the system with the FC holds over 
a quite wide region of energy up to £ ~ 6f- This means that beyond the FC phase transition 
point there is no room for the Landau theory in the vicinity of the origin £ = 0. Otherwise the 
ratio rc{e 0) should vanish and then some completely fiat portion ^ = in the spectrum 
of the single particle excitations should survive. But in this case calculations along the same 
lines as done above results in the nonzero value of rc{e — > 0). This contradiction proves the 
inapplicability of the Landau theory. 

Now we return to the damping 7„. From eq. ( ^1]) one has 



l2n{e) = C,r. I j^^^de, (69) 
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with the constant C2n whose value is less than 1. In principle, this integral is a linear in e 
function in the region e ^ 6p so that the total damping 7„(e) is 

Inie) = -ln£, £ < Sf, (70) 

the value of the constant 7„, proportional to ^^-20/^22 being, probably, considerably less than 
1. 

To establish the structure of the spectrum of the single particle excitations beyond the 
FC phase transition point one needs to solve self-consistently the system (^) and (^) 
accounting for the energy dependence of the damping 7(e) given by eqs. (|66D, (|6^) and ([T0|). 
One of the ways to obtain the solution is the iterative procedure. We start with the FC 
distribution n{p) found from eq. (|25| ) assuming the explicit form of this equation, where the 
damping is completely ignored, to be known. With this distribution at hand, we evaluate 
the parameters determining the damping at different energies. Having found these constants 
we substitute them into eq. (p8|), evaluate the new distribution function n(p) and insert it 
into eq. ( PT| ) to find the new single particle spectrum and then repeat all the procedure till 
its convergence. 

It should be kept in mind that the results depend on rj as an input parameter. What 
actually happens is that the FC density Pc is not the input parameter. It is given by an 
integral 

Pc = J n{p) dr. (71) 
n 

The integration in (|7lD is performed over regions where the slope of the single particle spec- 
trum is small compared to the usual one. In the quasiparticle approximation, the boundaries 
of this region are determined unambiguously, since this slope is identically 0, and the density 
Pc is also evaluated unambiguously being proportional to the excess — Tcr- However, the 
damping implicitly influences on the size of the FC region so that the integral (|7T1) turns 
out to be an equation for obtaining p^ and with it the set of equations describing the basic 
properties of the system beyond the FC phase transition point becomes complete. 

The numerical solution of this system is beyond of the scope of this article. We restrict 
ourselves to several comments. Since the ratio 'yj^e)!^ exceeds 1 from to bp the kink at 
,^ = in the distribution evaluated in the quasiparticle approximation is split: its upper 
part moves to the left while the low one, to the right and the magnitude of this shift is of 
order of bp. Steep pieces of the curve n(^) are connected by a bridge the length of which 
and the slope depend on input parameters. If it is almost horizontal, then the insertion of 
such a distribution n(^) into eq. (|25|) results in splitting the plateau in the single particle 
spectrum, one part of which is shifted upward while the other, downward of the Fermi level 
resembling a zigzag. The lower part occupied by electrons can be investigated in ARPES 
studies while the upper one, not. In this case, the spectrum of the former FC states remains 
practically fiat but it is shifted from the Fermi surface. Its value ~ bp can be treated as the 
characteristic FC energy Spc- 

On the other hand, if the slope of the bridge is not so small then the zigzag gets smeared. 
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Then the average slope of the spectrum is estimated from an approximate formula 



fd^. di . 
Sfc ^ —dpc::^—L. 
J dp dp 

n 

In the two-dimensional problem, the characteristic size L of the FC region is proportional 
to 772, in the three-dimensional system it is proportional to r]^. Substitution here epc — 
yields in the two-dimensional case the result d^/dp ~ 772 while in the three-dimensional 

2 

system d^/dp r]3. Thus we infer that the plateau in the single particle spectrum is inclined 
but the slope turns out to be small: it vanishes together with the condensate density pc- 



5. Comparison with the phenomenological theory by Varma et al. [pT[| 

Now we are able to compare the results obtained with postulates of the phenomenological 
theory by Varma et al. |^ of normal states of high-T^ superconductors. The first postulate 



is concerned with the polarizability P(q, a;) and requires it to be an cu-independent quantity 
in contrast to the ordinary linear tu-dependence. The second is related to the mass operator 



S(p,e). It is assumed to have the form (27) inherent in marginal Fermi liquid. As we have 
seen, such a behavior of the mass operator S beyond the FC phase transition point is due 
to flattening the single particle spectra. 

As for the polarizability P(q, a;) which determines a relation between an external field 
Vo and a real one V acting inside of the system: \^(q, cj) = P(q, a;) Vo(q, cj), we evaluate 
it neglecting connections between the normal and the FC regions, i.e. suggesting = 0. 
Then Jig P"Hq,o^) = 1 - J^(q)A„(q, oj) and 

ImP(g,o;) = ImA„(g,^) \P{qM? ^ Ar„(0) .^^^f,^ (72) 

In the optical conductivity problem, the value of the wave vector g = cj/c is so small that 
the damping effects exhibit themselves in full force. In the normal region, 

1 fj 'jn{£hn{£-^)dedp^dS 



(73) 

This integral is evaluated analytically and the result 

/ds 
I . ^|^| I ^- (74) 

is drastically enhanced compared to the conventional one (|4l|) , proportional to uj. Upon 
substituting into (^) the result 7„(£:) = — 7n£^ and simple integration one finds 

hn{q^Q.uj)--— (75) 

In 
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where 7„ ~ 1. With this result, simple algebra leads to the following estimate 

ImP(g<g,„a;)~-iVo(0). (76) 

Thus the first postulate of the theory |2T|] consisting of the independence of the polarizability 
P{uj) of the frequency u is reproduced in our analysis, too. 

With increasing u the damping grows up and the magnitude of ImP drops. As u 
increases from epc the damping 7(5) rapidly drops. As a result, the value of 6„(g — >• 0,a;) 
grows up. However, simultaneously the denominator of ([72|) increases, even faster than the 
numerator, and ImP(g 0,uj) falls rapidly to values of order of ?7A'"o(0). As a result, the 
optical conductivity of the system with the FC has the maximum at frequencies uj ~ rje^p. 
Thus in the systems with the FC there exists a new mechanism, different from the polaron 
one, for the appearance of irregular behavior of the optical conductivity a^u). 

The behavior of the real part a„(g, u) in the static limit also differs from the Fermi liquid 
one. Indeed 

a„(g=0,a;=0) ~ — -— / / (77) 

In contrast to the quasiparticle-approximation formula, this integrand has no singularities 
and therefore replacing integration over p by that over t = e— (T„(e)— ^„(p) and assuming the 
derivative dt/dp = —d^n{p)/dp to be constant one finds 

ReA{q=0,uj=0) =0. (78) 

The Landau result is restored only if g > g^- Of course, this result holds only if the 
momentum dependence of the group velocity ^„(p) is negligible but usually this dependence 
is sufficiently weak that results in a significant suppression of the real part of the particle-hole 
propagator and, in its turn, in the suppression of the well known results for the spin and 
charge susceptibility. From these results one obtains that the well known singular behavior 
of the propagator A{q,uj) in the static limit g — 0, — is absent in the systems with 
the FC and the reason for that is the linear in e damping of the single particle degrees of 
freedom. 



6. Conclusion 

In this article, proceeding from the standard many-body-theory formalism we analyzed 
the interplay between the fiattening and the damping in the single particle spectra of the 
strongly correlated Fermi systems beyond the FC phase transition point assuming the FC 
density pc to be relatively small: rj = pd p < 1. We concentrated on crude features of 
this phenomenon not caring about numerical factors in quantities of our interest. We have 
found that in the quasiparticle approximation there exists a class of solutions of the many 
body theory equation for the spectrum ^(p) with the zero value of the quasiparticle group 
velocity at the Fermi surface. The necessary condition for the appearance of these solutions 
is a sufficiently strong momentum dependence of the mass operator S(p,£:). Such a strong 
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momentum dependence may stem from the fact that in the vicinity of many phase transi- 
tion points, e.g. antiferromagnetism the effective interaction r(q) has quite narrow peaks. 
Employing the RPA for the evaluation of the contribution of these resonance components 
into the mass operator S(p, e) a well pronounced momentum dependence does emerge in 
this function. Thus we can infer that in strongly correlated electron systems of solids, the 
fermion condensation could be entailed by the spin fluctuations and, hence, the fermion 
condensation is the precursor of antiferromagnetism. 

We demonstrated that in the vicinity of the Fermi surface there is no room for the Landau 
theory: the system with the FC does behave itself as marginal Fermi liquid, i.e. the damping 
7(5) of the states adjacent to the Fermi surface is linear in energy. The enhancement of 7 
should entail a significant spread of the FC states and influences on the intensity of the peaks 
observed in ARPES studies in the region occupied by the FC. This intensity is proportional 

to the factor (^1— '^^g^'^-' l£=£(p)) therefore it is suppressed in the FC region adjacent to 

the van Hove points. Of course, the magnitude of this suppression greatly depends on input 
parameters and varies from one compound to another. 

Two different marginal regimes of the energy dependence of 7 have been established. The 
perturbative regime 7(e) ~ rjs is adequate at energies e > dp — rjEp. The strong coupling 
regime with 7(e) ~ e is involved in the vicinity of e = 0. The transition from the first to the 
second one is, probably, quite smooth, but dealing with the states near the van Hove points 
this transition could be well pronounced providing an unaccustomed shape of the energy 
dependence of 7(e). 

We established that behavior of the imaginary part of the polarization operator 
P{q = 0,uj) of the system with the FC drastically differs from that assigned by the Fermi 
liquid theory: its magnitude remains constant over the low-frequency region. This effect 
can lead to anomalies in the optical conductivity cr(to') conventionally associated with the 
polaron mechanism. We also uncovered an alteration of the real part of the particle-hole 
propagator that entails a suppression of the spin and charge susceptibility of the strongly 
correlated systems. 

It is worth noting that the presence of the marginal term in the mass operator should 
change many formulas of the BCS theory of superfiuid systems with the FC including the 
BCS relation between A(T = 0) and the critical temperature T^. This problem will be 
investigated in the next article in more detail. 

In conclusion, that the analysis of many experimental data on the base of Luttinger 
liquid where the renormalization constant z = is an appropriate guide to understanding 
of behavior of stronglycorrelated electron systems became clear, probably, ten years ago 
2l| , p3| . At the same time, many properties of these systems are known to be succesfully 



treated proceeding from Fermi liquid theory. What we demonstated in our article is that 
the fermion condensation model serves as a "bridge" connecting both these approaches. 
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Figures captions 



Fig. 1. Graphical solution of eq. ([T^). The distribution n{^) (bold line) is crossed by the 
set of straight lines — ^p) /V . The crossing points of these lines with the vertical line ^ = 
are included into the FC region Q. 

Fig. 2. The graphical representation of the momentum dependent contribution to the 
mass operator ([l9|) resulting from exchange spin fluctuation term (plSf ) of the effective inter- 
action. 

Fig. 3. FC regions in the square lattice within the model with the spin-spin effective 
interaction ( p!8|) peaked at the vector Q = {71,11). The formfactors 0(p) are assumed to 
vanish far from the van Hove points. The parameters used are: t=0.7, ^=l.leot, A=0.4£:o- 

Fig. 4. The residue of the single-particle Green function in marginal Fermi liquid with 
K{e) = — e as a function of energy e measured in 6%. 

Fig. 5. The graphical representation for the contribution 71^ to the total damping. 

Fig. 6. The same for 71c (a), 72„ (b) and 72c (c). 
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